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Abstract 

We have derived the integral form of the cusp and node coalescence conditions satisfied by the 
wavefunction at the coalescence of two charged particles in D > 2 dimension space. From it we 
have obtained the differential form of the coalescence conditions. These expressions reduce to 
the well-known integral and differential coalescence conditions in D = 3 space. It follows from 
the results derived that the approximate Laughlin wavefunction for the fractional Quantum Hall 
Effect satisfies the node coalescence condition. It is further noted that the integral form makes 
evident that unlike the electron- nucleus coalescence condition, the differential form of the electron- 
electron coalescence condition cannot be expressed in terms of the electron density at the point 
of coalescence. From the integral form, the integral and differential coalescence conditions for the 
pair-correlation function in D > 2 dimension space are also derived. The known differential form 
of the pair function cusp condition for the uniform electron gas in dimensions D = 2, 3 constitute 
a special case of the result derived. 
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I. INTRODUCTION 



The non relativistic wavefunction exhibits a cusp or a node at the coalescence of any 
two charged particles because of the singularity associated with the Coulomb interaction 
between them. The differential form of the cusp coalescence condition satisfied by the 
wavefunction for both electron-electron and electron-nucleus coalescence was derived for 
dimensions D = 3 by Kato yj. This form of the cusp condition is in terms of the spherical 
average of the wavefunction taken about the singularity. It is this differential form that is 
most extensively employed in the literature. The integral form of the electron-electron and 
electron-nucleus cusp and node coalescence condition for D = 3 was originally conjectured 
by BingelQl. They have been derived rigorously from the Schrodinger equation by Pack and 
Byers Brown |3(. The integral form of the coalescence conditions, however, is more general, 
and therefore more significant, because it retains the angular dependence of the wavefunc- 
tion at coalescence. The differential form is then readily obtained from the integral form 
by spherically averaging and differentiating. An example of the usefulness of the integral 
brm is the following. Over the past two and a half decades, there has been uncertainty 
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about whether the local effective potential energy of Kohn-Sham density functional 
theory was finite or singular at the nucleus of atoms. Employing the differential form 
of the cusp condition, it was initially proved |4j that the potential energy is finite at the 
nucleus of spherically symmetric atoms. However, using the integral form, the conclusion 
of finiteness can be generalized p to atoms of arbitrary symmetry, to molecules, and to 
periodic and aperiodic solids. As another example, it is well known that the differential 
form of the electron-nucleus cusp condition may be expressed j?| in terms of the spherical 
average of the electron density about the nucleus. It has consequently been assumed Q] that 
the differential form of the electron-electron cusp condition can similarly be expressed in 
terms of the density at the point of electron coalescence. The integral form of the electron- 
electron cusp condition makes clear that no such expression in terms of the density can exist. 

At present, there is considerable interest in both lower dimensional systems such as 
the two-dimensional (D = 2) electron gas [l| , the two-dimensional gas in the presence of 
a magnetic field ( the Quantum Hall Effect) [10(, and in their higher-dimensional D > 4 



generalizations llj. Motivated by this interest, we have derived the integral form of the 
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cusp and node coalescence condition for dimensions D > 2. It becomes evident that the 
approximate Laughlin[lo| wavefunction for the fractional Quantum Hall effect satisfies the 
D = 2 node coalescence condition. The D > 2 differential form of the cusp condition is 
then obtained by spherically averaging about the singularity over D-dimensional space and 
differentiating. Employing the integral cusp condition, we have also derived the integral 
and differential forms of the pair correlation function at the coalescence of two identical 
particles in D > 2 dimensions. The differential form of the cusp condition for the pair 
correlation function for the uniform electron gas for D = 3 3| and D = 2^| constitute a 
special case of our result. The use of cusp conditions for the pair-correlation function in 
D = 3 spin-density functional theory is expounded in Ref. 14. More recent work on nuclear 
cusp conditions relevant to density functional theory is described in Ref. 15. 



II . DERIVATION 

Consider a nonrelativistic system of N charged particles in D (> 2) dimension space with 
the Hamitonian (h = e = 1) 

N i N 7 7 

» = -Y,^+Y,^ (i) 

i=l *"H j>i=l ' % 3 

where wij and Zi are the mass and charge of the ith particle, and = |r^ — In D 



dimesion space, r = (xi, x 2 , ...id), r = yJ2k=i x k x k an d V 2 = J2k=i ^"af - - -^ ue t° the 
Coulomb potential energy term, the Hamiltonian is singular when any two particle i and 
j coalesce^- — > 0). For the wavefunction \&(ri, r2, ...r/v) which satisfies the Schrodinger 
equation 

r 2 , ...r N ) = E^(r u r 2 , ...r N ), (2) 

to be bounded and remain finite at the singularites, it must satisfy a cusp coalescence 
condition. If the wavefunction vanishes at the singularity, it must satisfy a node coalescence 
condition. 

We are interested in the form of the wavefunction when two particles approach each other, 
i.e. when is very small. Following Pack and Byers Brown 3] we focus our attention on 
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two particles 1 and 2, and transform their coordinates ri and r 2 to the center of mass Ri 2 
and relative coordinates r i2 as 

m 1 r 1 + m 2 r 2 

*M 2 — ; , 

mi + m 2 

r 12 = n - r 2 . (3) 
The Hamiltonian of Eq.(l) may then be rewritten as 

AT -I AT 7 7 

+E^v;, + £^, ( 4) 

i=3 zm * j>«=3 '*3 

where fii 2 = mi+m 2 * s tne re duced mass of particles 1 and 2. When particle 1 and 2 are 
within a small distance of each other ( < r 12 < e), and all other particles are well separated, 
there is only one singularity in the Hamiltonian. Retaining only terms of lower order in r 12 , 
Eq. (3) reduces to 

[-^-V2 U + ^ + 0(e°)]*(ri, r 2) ...r,) = 0, (5) 
^12 r 12 

where 0(e°) refers to terms of order zero (constant) and higher order in r i2 and the vector 
components (r 12 ) fc , k — 1, 2, ...D). 

Eq.(4) is a one-electron-atom equation in D dimension space. Furthermore, it is not 
an eigenvalue equation. We are interested in the solution of the equation that are finite 
and continuous at the singularities. Thus, in the limit as ri — > r 2 , we could write the 
wavefunction in two parts as 

(ri, r 2 , ...rjv) = *(r 2 , r 2 , r 3 ..., r^) + <f*(ri, r 2 , ...1^), (6) 

where the term <5\&(ri, r 2 , ...rjv) must vanish at the singularity ri = r 2 . From the above 
it follows that we need consider only terms of first-order in r 12 in 5^ when ri is near r 2 . 
Dropping the subscript 12, and writing r 12 as r, we see that 5^f could have terms of the 
form r D B(r 2 ,r 3 ...,r N ), or Y,i fiBi(r 2 , ...r N ) where I denotes the 2 to (D — l)-dimensional 
subspace and 77 the distance in the subspace, or r ■ C(r 2 , ...rjv). Here, td — r is just 
the conventional distance in D dimensional space. As an example of the second kind of 



terms, consider the case of D= 4. There are four 3-dimensional subspaces, and six 2- 
dimensional subspaces. The distances in the 3-dimensional subspaces constitute the terms 



\Jx\ + x\ + £3, \jx\ + x\ + x\,\Jx\ + x 3 + x|, and ^/xf + x 2 + x|. The distances in the 2- 



dimensional subspaces are yx 2 + xj, where i 7^ ,7, and i,j = 1,2,3,4. In the limit r — > 0, 
such terms do not affect the 0(D) symmetry of the first two terms on the left side of Eq.(5). 
However, the Laplacian of these terms is singular since V 2 r<j = where is the distance 
in the d-dimensional subspace. These singularities cannot be cancelled by the Coulomb 
potential energy singularity term of Eq.(5). Thus, their coefficients must vanish. Therefore, 
(recovering the subscript 12), Eq.(6) must be written as 

#(ri, r 2 , ...r;v) = *(r 2 , r 2 , r 3 ..., r N ) + r 12 B(r 2 , r 3 , r N ) + r 12 • C(r 2 , r 3 , r N ) + 0(e 2 ) (7) 

Next substitute Eq. (6) for ^(ri, r 2 , ...rjy) into Eq.(4) and employ V 2 r = 0. For the Coulomb 
potential energy singularity to be cancelled we have 

1 D-l , . Z X Z 2 , 

-B(r 2 , r 3 , rjv) + *(r 2 , r 2 , r 3 , rjv) = 0. (8) 



2//; 2 ri2 " ' " ~ r 12 
At the point of coalescence, the wavefunction is 4 , (r 2 , r 2 , r 3 , r^)- Since we require the 
wavefunction to be finite at this point we have 

B(r 2 , r 3 , r N ) = 2Z ^ Z ^ 12 ^(r 2 , r 2 , r 3 ..., r^). (9) 

Thus, in the limit as T\ 2 — > 0, we may write the wavefunction as 

*(ri,r 2 , ...r N ) = #(r 2 ,r 2 ,r 3 , ...,r N )(l + ^ _ 2 ^ 12 r i2 ) + r 12 • C(r 2 ,r 3 , ...,rjv). (10) 

This is the general form of the integral cusp coalesence condition in D-dimension space. Note 
that this expression is equally valid even if the wavefunction vanishes at the point of coales- 
cence, i.e. if \l/(r 2 ,r 2 , ...,rjy) = 0. This latter is referred to as a node coalescence condition. 
By taking the spherical average of Eq.(10) about the point of coalescence in D-dimension 
space, we obtain the differential form of the cusp condition: 

2Z 1 Z m2 

(^HIm-> 0=-^- *(ri2 = 0), (11) 

where ^ is the spherically averaged wavefunction. For the electron-nucleus coalescence, 
Z\ = —1,Z 2 = Z the nuclear charge , and n\ 2 ~ m e the mass of the electron. For the 
electron-electron coalescence, Z\ = —1,Z 2 = — l,/^i 2 = \m e . In D = 3 dimensions, the 



traditional integral and differential cusp conditions are recovered. 



It follows from the definition of the density p(r x ) which is 

p{r x ) = N J **(n, r 2 , ...r N )^(r u r 2 , ...r N )dr 2 ...dr N , (12) 

that the differential form of the electron-nucleus cusp condition may be expessed in terms 
of the density at and about the nucleus. In this case with the nuclear positions fixed, 
the wavefunction does not depend upon the nuclear coordinates. Then when the electron 
at ri approaches a nucleus, say fixed at the origin, the integral for the density which is 
over all the remaining N — 1 electrons can be performed. It is evident from the integral 
cusp coalesence condition Eq.(lO) and the above definition of the density, that there can 
be no such differential form of the cusp condition for electron- electron coalescence. Such 
expressions in the literature |8( are therefore invalid. 

It is however, possible to derive the integral and differential forms of the cusp condition 
in D-dimensions for the two-particle correlation function p(rir 2 ) defined as 

N(N-l) f 

p{r 1 r 2 ) = J #*(ri, r 2 , ...rjv)*(ri, r 2 , ...r N )dr 3l ...dr N . (13) 

In the limit as r x — > r 2 and employing the integral cusp expression of Eq.(lO) we have that 
to first-order in r± 2 

p(r x r 2 ) = p(r 2 r 2 )(l + 2Z ^ Z ^ U y 12 )2 + — — J (r 12 ■ C(r 2 , ...r N )) 2 dr 3 ...dr N 

+ N(N-l)(l+ 2Z ^ i U r 12 ) J(v 12 -C(v 2 ,...v N ))^(r 2 ,r 2 ,...r N )dr 3 ...dv N 

+ 0(e 3 ). (14) 

This is the integral form of the cusp coalescence condition for the pair- correlation function. 
The differential form of the cusp condition for the correlation function is 

<9p(r x r 2 ) AZtZ 2 fi 12 



v dr 12 Ari2 ~*° D-l 
where p is the spherical average of the correlation function. The differential cusp condition 
derived for the uniform electron gas in dimensions D = 3 and 2 constitute a special case of 
this general form of the differential cusp condition. 



p(r 2 r 2 ). (15) 
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III. CONCLUDING REMARKS 



In this paper we have derived the integral and differential forms of the coalescence con- 
dition satisfied by the wavefunction at the coalescence of two charged particles in D > 2 
dimensions. The corresponding integral and differential cusp conditions for the pair corre- 
lation function in D > 2 dimensional space are thereby derived. These results are also valid 
for the case of a more generalized Hamitonian which includes a term of the form J2iv(ri), 
where v (r*) is a local external potential. This is because such terms do not contribute any 
singularities at the coalescense of the two charged particles. Such an external potential could 
effectively arise for example in the presence of an electric or magnetic field. The coalescence 
conditions constitute rigorous constraints on the construction and evaluation of approximate 
wavef unctions. Thus, having derived the integral coalesence condition for D = 2, we now 
understand that the approximate Laughlin wavefunction does in fact satisfy the node 
coalescence condition. The Laughlin wavefunction is 

N i N 

* m (z u ...,z N ) = Hizj - z k ) m exp[--J2\zi\ 2 }, (16) 
j<k 4 i 

where m is an odd integer, and Zj = Xj + iyj is the location of the j th particle expressed 
as a complex number. It is evident that this wavefunction satisfies the node coalescence 
condition for m > 3. The cusp conditions also impose constraints on the construction 



0,0, 



and on the 



of approximate energy functionals within density functional theory 
wavefunction functionals of the density employed in quantal density functional theory 
The integral cusp condition in one dimension, (which is different in form from the result 
derived for dimension D > 2), and from wh.eh emanates some very interesting phys.es, is 
to be published elsewhere [181]. 

This work was supported in part by the Research Foundation of CUNY. 
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